3.601 \(\int \sqrt {\sec (c+d x)} (a+b \sec (c+d x))^4 \, dx\)

Optimal. Leaf size=247 \[ \frac {2 b^2 \left (39 a^2+5 b^2\right ) \sin (c+d x) \sec ^{\frac {3}{2}}(c+d x)}{21 d}+\frac {8 a b \left (5 a^2+3 b^2\right ) \sin (c+d x) \sqrt {\sec (c+d x)}}{5 d}-\frac {8 a b \left (5 a^2+3 b^2\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} E\left (\left .\frac {1}{2} (c+d x)\right |2\right )}{5 d}+\frac {2 \left (21 a^4+42 a^2 b^2+5 b^4\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} F\left (\left .\frac {1}{2} (c+d x)\right |2\right )}{21 d}+\frac {36 a b^3 \sin (c+d x) \sec ^{\frac {5}{2}}(c+d x)}{35 d}+\frac {2 b^2 \sin (c+d x) \sec ^{\frac {3}{2}}(c+d x) (a+b \sec (c+d x))^2}{7 d} \]

[Out]

2/21*b^2*(39*a^2+5*b^2)*sec(d*x+c)^(3/2)*sin(d*x+c)/d+36/35*a*b^3*sec(d*x+c)^(5/2)*sin(d*x+c)/d+2/7*b^2*sec(d*
x+c)^(3/2)*(a+b*sec(d*x+c))^2*sin(d*x+c)/d+8/5*a*b*(5*a^2+3*b^2)*sin(d*x+c)*sec(d*x+c)^(1/2)/d-8/5*a*b*(5*a^2+
3*b^2)*(cos(1/2*d*x+1/2*c)^2)^(1/2)/cos(1/2*d*x+1/2*c)*EllipticE(sin(1/2*d*x+1/2*c),2^(1/2))*cos(d*x+c)^(1/2)*
sec(d*x+c)^(1/2)/d+2/21*(21*a^4+42*a^2*b^2+5*b^4)*(cos(1/2*d*x+1/2*c)^2)^(1/2)/cos(1/2*d*x+1/2*c)*EllipticF(si
n(1/2*d*x+1/2*c),2^(1/2))*cos(d*x+c)^(1/2)*sec(d*x+c)^(1/2)/d

________________________________________________________________________________________

Rubi [A]  time = 0.37, antiderivative size = 247, normalized size of antiderivative = 1.00, number of steps used = 9, number of rules used = 8, integrand size = 23, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.348, Rules used = {3842, 4076, 4047, 3768, 3771, 2639, 4046, 2641} \[ \frac {2 b^2 \left (39 a^2+5 b^2\right ) \sin (c+d x) \sec ^{\frac {3}{2}}(c+d x)}{21 d}+\frac {8 a b \left (5 a^2+3 b^2\right ) \sin (c+d x) \sqrt {\sec (c+d x)}}{5 d}+\frac {2 \left (42 a^2 b^2+21 a^4+5 b^4\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} F\left (\left .\frac {1}{2} (c+d x)\right |2\right )}{21 d}-\frac {8 a b \left (5 a^2+3 b^2\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} E\left (\left .\frac {1}{2} (c+d x)\right |2\right )}{5 d}+\frac {36 a b^3 \sin (c+d x) \sec ^{\frac {5}{2}}(c+d x)}{35 d}+\frac {2 b^2 \sin (c+d x) \sec ^{\frac {3}{2}}(c+d x) (a+b \sec (c+d x))^2}{7 d} \]

Antiderivative was successfully verified.

[In]

Int[Sqrt[Sec[c + d*x]]*(a + b*Sec[c + d*x])^4,x]

[Out]

(-8*a*b*(5*a^2 + 3*b^2)*Sqrt[Cos[c + d*x]]*EllipticE[(c + d*x)/2, 2]*Sqrt[Sec[c + d*x]])/(5*d) + (2*(21*a^4 +
42*a^2*b^2 + 5*b^4)*Sqrt[Cos[c + d*x]]*EllipticF[(c + d*x)/2, 2]*Sqrt[Sec[c + d*x]])/(21*d) + (8*a*b*(5*a^2 +
3*b^2)*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/(5*d) + (2*b^2*(39*a^2 + 5*b^2)*Sec[c + d*x]^(3/2)*Sin[c + d*x])/(21*d
) + (36*a*b^3*Sec[c + d*x]^(5/2)*Sin[c + d*x])/(35*d) + (2*b^2*Sec[c + d*x]^(3/2)*(a + b*Sec[c + d*x])^2*Sin[c
+ d*x])/(7*d)

Rule 2639

Int[Sqrt[sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2*EllipticE[(1*(c - Pi/2 + d*x))/2, 2])/d, x] /; FreeQ[{
c, d}, x]

Rule 2641

Int[1/Sqrt[sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[(2*EllipticF[(1*(c - Pi/2 + d*x))/2, 2])/d, x] /; FreeQ
[{c, d}, x]

Rule 3768

Int[(csc[(c_.) + (d_.)*(x_)]*(b_.))^(n_), x_Symbol] :> -Simp[(b*Cos[c + d*x]*(b*Csc[c + d*x])^(n - 1))/(d*(n -
1)), x] + Dist[(b^2*(n - 2))/(n - 1), Int[(b*Csc[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1
] && IntegerQ[2*n]

Rule 3771

Int[(csc[(c_.) + (d_.)*(x_)]*(b_.))^(n_), x_Symbol] :> Dist[(b*Csc[c + d*x])^n*Sin[c + d*x]^n, Int[1/Sin[c + d
*x]^n, x], x] /; FreeQ[{b, c, d}, x] && EqQ[n^2, 1/4]

Rule 3842

Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_), x_Symbol] :> -Simp[(b^2*
Cot[e + f*x]*(a + b*Csc[e + f*x])^(m - 2)*(d*Csc[e + f*x])^n)/(f*(m + n - 1)), x] + Dist[1/(d*(m + n - 1)), In
t[(a + b*Csc[e + f*x])^(m - 3)*(d*Csc[e + f*x])^n*Simp[a^3*d*(m + n - 1) + a*b^2*d*n + b*(b^2*d*(m + n - 2) +
3*a^2*d*(m + n - 1))*Csc[e + f*x] + a*b^2*d*(3*m + 2*n - 4)*Csc[e + f*x]^2, x], x], x] /; FreeQ[{a, b, d, e, f
, n}, x] && NeQ[a^2 - b^2, 0] && GtQ[m, 2] && (IntegerQ[m] || IntegersQ[2*m, 2*n]) &&  !(IGtQ[n, 2] &&  !Integ
erQ[m])

Rule 4046

Int[(csc[(e_.) + (f_.)*(x_)]*(b_.))^(m_.)*(csc[(e_.) + (f_.)*(x_)]^2*(C_.) + (A_)), x_Symbol] :> -Simp[(C*Cot[
e + f*x]*(b*Csc[e + f*x])^m)/(f*(m + 1)), x] + Dist[(C*m + A*(m + 1))/(m + 1), Int[(b*Csc[e + f*x])^m, x], x]
/; FreeQ[{b, e, f, A, C, m}, x] && NeQ[C*m + A*(m + 1), 0] &&  !LeQ[m, -1]

Rule 4047

Int[(csc[(e_.) + (f_.)*(x_)]*(b_.))^(m_.)*((A_.) + csc[(e_.) + (f_.)*(x_)]*(B_.) + csc[(e_.) + (f_.)*(x_)]^2*(
C_.)), x_Symbol] :> Dist[B/b, Int[(b*Csc[e + f*x])^(m + 1), x], x] + Int[(b*Csc[e + f*x])^m*(A + C*Csc[e + f*x
]^2), x] /; FreeQ[{b, e, f, A, B, C, m}, x]

Rule 4076

Int[((A_.) + csc[(e_.) + (f_.)*(x_)]*(B_.) + csc[(e_.) + (f_.)*(x_)]^2*(C_.))*(csc[(e_.) + (f_.)*(x_)]*(d_.))^
(n_.)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)), x_Symbol] :> -Simp[(b*C*Csc[e + f*x]*Cot[e + f*x]*(d*Csc[e + f*x
])^n)/(f*(n + 2)), x] + Dist[1/(n + 2), Int[(d*Csc[e + f*x])^n*Simp[A*a*(n + 2) + (B*a*(n + 2) + b*(C*(n + 1)
+ A*(n + 2)))*Csc[e + f*x] + (a*C + B*b)*(n + 2)*Csc[e + f*x]^2, x], x], x] /; FreeQ[{a, b, d, e, f, A, B, C,
n}, x] &&  !LtQ[n, -1]

Rubi steps

\begin {align*} \int \sqrt {\sec (c+d x)} (a+b \sec (c+d x))^4 \, dx &=\frac {2 b^2 \sec ^{\frac {3}{2}}(c+d x) (a+b \sec (c+d x))^2 \sin (c+d x)}{7 d}+\frac {2}{7} \int \sqrt {\sec (c+d x)} (a+b \sec (c+d x)) \left (\frac {1}{2} a \left (7 a^2+b^2\right )+\frac {1}{2} b \left (21 a^2+5 b^2\right ) \sec (c+d x)+9 a b^2 \sec ^2(c+d x)\right ) \, dx\\ &=\frac {36 a b^3 \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d}+\frac {2 b^2 \sec ^{\frac {3}{2}}(c+d x) (a+b \sec (c+d x))^2 \sin (c+d x)}{7 d}+\frac {4}{35} \int \sqrt {\sec (c+d x)} \left (\frac {5}{4} a^2 \left (7 a^2+b^2\right )+7 a b \left (5 a^2+3 b^2\right ) \sec (c+d x)+\frac {5}{4} b^2 \left (39 a^2+5 b^2\right ) \sec ^2(c+d x)\right ) \, dx\\ &=\frac {36 a b^3 \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d}+\frac {2 b^2 \sec ^{\frac {3}{2}}(c+d x) (a+b \sec (c+d x))^2 \sin (c+d x)}{7 d}+\frac {4}{35} \int \sqrt {\sec (c+d x)} \left (\frac {5}{4} a^2 \left (7 a^2+b^2\right )+\frac {5}{4} b^2 \left (39 a^2+5 b^2\right ) \sec ^2(c+d x)\right ) \, dx+\frac {1}{5} \left (4 a b \left (5 a^2+3 b^2\right )\right ) \int \sec ^{\frac {3}{2}}(c+d x) \, dx\\ &=\frac {8 a b \left (5 a^2+3 b^2\right ) \sqrt {\sec (c+d x)} \sin (c+d x)}{5 d}+\frac {2 b^2 \left (39 a^2+5 b^2\right ) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{21 d}+\frac {36 a b^3 \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d}+\frac {2 b^2 \sec ^{\frac {3}{2}}(c+d x) (a+b \sec (c+d x))^2 \sin (c+d x)}{7 d}-\frac {1}{5} \left (4 a b \left (5 a^2+3 b^2\right )\right ) \int \frac {1}{\sqrt {\sec (c+d x)}} \, dx+\frac {1}{21} \left (21 a^4+42 a^2 b^2+5 b^4\right ) \int \sqrt {\sec (c+d x)} \, dx\\ &=\frac {8 a b \left (5 a^2+3 b^2\right ) \sqrt {\sec (c+d x)} \sin (c+d x)}{5 d}+\frac {2 b^2 \left (39 a^2+5 b^2\right ) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{21 d}+\frac {36 a b^3 \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d}+\frac {2 b^2 \sec ^{\frac {3}{2}}(c+d x) (a+b \sec (c+d x))^2 \sin (c+d x)}{7 d}-\frac {1}{5} \left (4 a b \left (5 a^2+3 b^2\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \sqrt {\cos (c+d x)} \, dx+\frac {1}{21} \left (\left (21 a^4+42 a^2 b^2+5 b^4\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {1}{\sqrt {\cos (c+d x)}} \, dx\\ &=-\frac {8 a b \left (5 a^2+3 b^2\right ) \sqrt {\cos (c+d x)} E\left (\left .\frac {1}{2} (c+d x)\right |2\right ) \sqrt {\sec (c+d x)}}{5 d}+\frac {2 \left (21 a^4+42 a^2 b^2+5 b^4\right ) \sqrt {\cos (c+d x)} F\left (\left .\frac {1}{2} (c+d x)\right |2\right ) \sqrt {\sec (c+d x)}}{21 d}+\frac {8 a b \left (5 a^2+3 b^2\right ) \sqrt {\sec (c+d x)} \sin (c+d x)}{5 d}+\frac {2 b^2 \left (39 a^2+5 b^2\right ) \sec ^{\frac {3}{2}}(c+d x) \sin (c+d x)}{21 d}+\frac {36 a b^3 \sec ^{\frac {5}{2}}(c+d x) \sin (c+d x)}{35 d}+\frac {2 b^2 \sec ^{\frac {3}{2}}(c+d x) (a+b \sec (c+d x))^2 \sin (c+d x)}{7 d}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]  time = 1.72, size = 168, normalized size = 0.68 \[ \frac {2 \sec ^{\frac {7}{2}}(c+d x) \left (-84 a b \left (5 a^2+3 b^2\right ) \cos ^{\frac {7}{2}}(c+d x) E\left (\left .\frac {1}{2} (c+d x)\right |2\right )+b \sin (c+d x) \left (84 a \left (5 a^2+3 b^2\right ) \cos ^3(c+d x)+5 b \left (42 a^2+5 b^2\right ) \cos ^2(c+d x)+15 b^3\right )+5 \left (21 a^4+42 a^2 b^2+5 b^4\right ) \cos ^{\frac {7}{2}}(c+d x) F\left (\left .\frac {1}{2} (c+d x)\right |2\right )+42 a b^3 \sin (2 (c+d x))\right )}{105 d} \]

Antiderivative was successfully verified.

[In]

Integrate[Sqrt[Sec[c + d*x]]*(a + b*Sec[c + d*x])^4,x]

[Out]

(2*Sec[c + d*x]^(7/2)*(-84*a*b*(5*a^2 + 3*b^2)*Cos[c + d*x]^(7/2)*EllipticE[(c + d*x)/2, 2] + 5*(21*a^4 + 42*a
^2*b^2 + 5*b^4)*Cos[c + d*x]^(7/2)*EllipticF[(c + d*x)/2, 2] + b*(15*b^3 + 5*b*(42*a^2 + 5*b^2)*Cos[c + d*x]^2
+ 84*a*(5*a^2 + 3*b^2)*Cos[c + d*x]^3)*Sin[c + d*x] + 42*a*b^3*Sin[2*(c + d*x)]))/(105*d)

________________________________________________________________________________________

fricas [F]  time = 0.52, size = 0, normalized size = 0.00 \[ {\rm integral}\left ({\left (b^{4} \sec \left (d x + c\right )^{4} + 4 \, a b^{3} \sec \left (d x + c\right )^{3} + 6 \, a^{2} b^{2} \sec \left (d x + c\right )^{2} + 4 \, a^{3} b \sec \left (d x + c\right ) + a^{4}\right )} \sqrt {\sec \left (d x + c\right )}, x\right ) \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sec(d*x+c)^(1/2)*(a+b*sec(d*x+c))^4,x, algorithm="fricas")

[Out]

integral((b^4*sec(d*x + c)^4 + 4*a*b^3*sec(d*x + c)^3 + 6*a^2*b^2*sec(d*x + c)^2 + 4*a^3*b*sec(d*x + c) + a^4)
*sqrt(sec(d*x + c)), x)

________________________________________________________________________________________

giac [F]  time = 0.00, size = 0, normalized size = 0.00 \[ \int {\left (b \sec \left (d x + c\right ) + a\right )}^{4} \sqrt {\sec \left (d x + c\right )}\,{d x} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sec(d*x+c)^(1/2)*(a+b*sec(d*x+c))^4,x, algorithm="giac")

[Out]

integrate((b*sec(d*x + c) + a)^4*sqrt(sec(d*x + c)), x)

________________________________________________________________________________________

maple [B]  time = 12.66, size = 925, normalized size = 3.74 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(sec(d*x+c)^(1/2)*(a+b*sec(d*x+c))^4,x)

[Out]

-(-(-2*cos(1/2*d*x+1/2*c)^2+1)*sin(1/2*d*x+1/2*c)^2)^(1/2)*(2*a^4*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(-2*cos(1/2*d*x
+1/2*c)^2+1)^(1/2)/(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)*EllipticF(cos(1/2*d*x+1/2*c),2^(1/2))+
2*b^4*(-1/56*cos(1/2*d*x+1/2*c)*(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)/(-1/2+cos(1/2*d*x+1/2*c)^
2)^4-5/42*cos(1/2*d*x+1/2*c)*(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)/(-1/2+cos(1/2*d*x+1/2*c)^2)^
2+5/21*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(-2*cos(1/2*d*x+1/2*c)^2+1)^(1/2)/(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2
*c)^2)^(1/2)*EllipticF(cos(1/2*d*x+1/2*c),2^(1/2)))-8/5*a*b^3/(8*sin(1/2*d*x+1/2*c)^6-12*sin(1/2*d*x+1/2*c)^4+
6*sin(1/2*d*x+1/2*c)^2-1)/sin(1/2*d*x+1/2*c)^2*(12*EllipticE(cos(1/2*d*x+1/2*c),2^(1/2))*(2*sin(1/2*d*x+1/2*c)
^2-1)^(1/2)*(sin(1/2*d*x+1/2*c)^2)^(1/2)*sin(1/2*d*x+1/2*c)^4-24*cos(1/2*d*x+1/2*c)*sin(1/2*d*x+1/2*c)^6-12*El
lipticE(cos(1/2*d*x+1/2*c),2^(1/2))*(2*sin(1/2*d*x+1/2*c)^2-1)^(1/2)*(sin(1/2*d*x+1/2*c)^2)^(1/2)*sin(1/2*d*x+
1/2*c)^2+24*sin(1/2*d*x+1/2*c)^4*cos(1/2*d*x+1/2*c)+3*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(2*sin(1/2*d*x+1/2*c)^2-1)^
(1/2)*EllipticE(cos(1/2*d*x+1/2*c),2^(1/2))-8*sin(1/2*d*x+1/2*c)^2*cos(1/2*d*x+1/2*c))*(-2*sin(1/2*d*x+1/2*c)^
4+sin(1/2*d*x+1/2*c)^2)^(1/2)+8*a^3*b*(-(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)*(sin(1/2*d*x+1/2*
c)^2)^(1/2)*(2*sin(1/2*d*x+1/2*c)^2-1)^(1/2)*EllipticE(cos(1/2*d*x+1/2*c),2^(1/2))+2*(-2*sin(1/2*d*x+1/2*c)^4+
sin(1/2*d*x+1/2*c)^2)^(1/2)*cos(1/2*d*x+1/2*c)*sin(1/2*d*x+1/2*c)^2)/sin(1/2*d*x+1/2*c)^2/(2*sin(1/2*d*x+1/2*c
)^2-1)+12*a^2*b^2*(-1/6*cos(1/2*d*x+1/2*c)*(-2*sin(1/2*d*x+1/2*c)^4+sin(1/2*d*x+1/2*c)^2)^(1/2)/(-1/2+cos(1/2*
d*x+1/2*c)^2)^2+1/3*(sin(1/2*d*x+1/2*c)^2)^(1/2)*(-2*cos(1/2*d*x+1/2*c)^2+1)^(1/2)/(-2*sin(1/2*d*x+1/2*c)^4+si
n(1/2*d*x+1/2*c)^2)^(1/2)*EllipticF(cos(1/2*d*x+1/2*c),2^(1/2))))/sin(1/2*d*x+1/2*c)/(2*cos(1/2*d*x+1/2*c)^2-1
)^(1/2)/d

________________________________________________________________________________________

maxima [F(-1)]  time = 0.00, size = 0, normalized size = 0.00 \[ \text {Timed out} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sec(d*x+c)^(1/2)*(a+b*sec(d*x+c))^4,x, algorithm="maxima")

[Out]

Timed out

________________________________________________________________________________________

mupad [F]  time = 0.00, size = -1, normalized size = -0.00 \[ \int {\left (a+\frac {b}{\cos \left (c+d\,x\right )}\right )}^4\,\sqrt {\frac {1}{\cos \left (c+d\,x\right )}} \,d x \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a + b/cos(c + d*x))^4*(1/cos(c + d*x))^(1/2),x)

[Out]

int((a + b/cos(c + d*x))^4*(1/cos(c + d*x))^(1/2), x)

________________________________________________________________________________________

sympy [F(-1)]  time = 0.00, size = 0, normalized size = 0.00 \[ \text {Timed out} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(sec(d*x+c)**(1/2)*(a+b*sec(d*x+c))**4,x)

[Out]

Timed out

________________________________________________________________________________________